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Abstract 

We have investigated the influence of potassium intercalation on the formation of the charge- 
density wave (CD W) instability in 2i/-TaSe2 by means of Electron Energy-Loss Spectroscopy and 
density functional theory. Our observations are consistent with a filling of the conduction band as 
indicated by a substantial decrease of the plasma frequency in experiment and theory. In addition, 
elastic scattering clearly points to a destruction of the CDW upon intercalation as can be seen by 
a vanishing of the corresponding superstructures. This is accompanied by a new superstructure, 
which can be attributed to the intercalated potassium. Based on the behavior of the c-axis upon 
intercalation we argue in favor of interlayer-sites for the alkali-metal and that the lattice remains 
in the 2H modification. 

PACS numbers: 71.45.Lr,71.20.Tx,71.15.Mb 
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INTRODUCTION 



Due to their layered structure, transition-metal dichalcogenides (TMDCs) show a variety 
of anisotropic properties, such as thermal expansion, sound velocity and thermal as well as 
electrical conductivity [IHI]. The lowered dimensionality often also manifests itself in the 
occurrence of a charge-density wave (CDW) phase transition within the afe-plane of the crys- 
tal. In the present example of 2if-TaSe2, this CDW superstructure is incommensurate below 
120 K and becomes commensurate going below 90 K [5H7]. A detailed study of the phase 
transition — containing an incommensurate superstructure followed by a lock-in transition to 
a commensurate ordering vector — has for example been realised by Neutron scattering [7j 
or x-ray diffraction [8j. This transition as well as the CDW state itself were also investi- 
gated by angle- resolved photoemission measurements [9j [10] , clearly showing the opening of 
a gap at the Fermi level when going to the ordered phase. However, the reasoning that the 
Fermi surface is partially gapped already for temperatures far above the phase transition 
and the link of this fact to a possible nesting mechanism for the CDW order is supported 
by reflectivity as well as resistivity measurements [11J . 

Recent spectroscopic investigations on 2/7-TaSe2 showed results seemingly contradicting 
the conventional description of a metal — apart from the observation of the partial gapping — 
namely, a negative slope of the plasmon energy dispersion, which can be assigned to an 
interplay of the charge density fluctuations and the plasma resonance [12] . 

While the afr-planes of the crystals of the TMDCs consist of hexagonal layers, the c- 
direction is formed by sandwich-like structures of these layers, bound to each other through 
van-der- VKaa/s-interaction. The gaps between the building blocks are the typical place for 
all kinds of intercalates to arrange within the crystal [131 EH]- Thereby, the range of possible 
doping substances is as wide as the spectrum of physical properties to tune with intercalation. 

From the electronic point of view it is reported that the band structure changes extensively 
upon alkali metal (Na with Na/Ta= 0.3 and Cs with Cs/Ta varying between 0.3 and 0.6) 
doping in 2if-TaSe2 [IHl [16] and related compounds, such as lT-Rb x TaSe2 [H] and 1T- 
Cs x TaSe2 ( up to x = 0.5) [18J. 

Structural changes upon intercalation can be observed in various polytypes of the 
TMDCs. Primarily the c-axis is widened, while the spacing within the afr-plane remains 
nearly unchanged [HI [20] . Slight changes in the in-plane crystal structure indicated by a 
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change of super structure are observed for different compounds such as the lT-Cs x TaSe2 
mentioned before [IS] as well as in K x TiS 2 [2"T] . Because of the suppression of the CDW, 
as for example in 2H-K X T&S 2 [22], the observed super structures can be assumed to result 
from the intercalated atoms. 

Structural changes are also discussed in terms of staging effects in the class of silver doped 
TMDCs in close analogy to intercalated graphite compounds. The added silver atoms are 
supposed to slightly disarrange the chalcogenide atoms and with it opening new scattering 
channels. The observed stages are at doping rates of one and two thirds, respectively, for 
Aga;TaS2 as well as Ag x TiS2 |23j. Other examples also show super structure changes in 
Ag x NbS 2 [211125]. 

We report on the in situ intercalation into the bulk sample of 2if-TaSe2 with the alkali 
metal potassium, reaching a comparatively high stable doping rate. Electron Energy-Loss 
Spectroscopy is applied to investigate changes in the electronic as well as crystal structure. 
We demonstrate that doping levels up to K .77TaSe 2 can be reached which are characterized 
by well defined plasmon excitations, while the CDW superstructure is suppressed upon 
doping. 



EXPERIMENT AND CALCULATIONS 
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Figure 1. Charge density wave superstructure reflections arising by cooling down below the phase 
transition at 120 K. The figures show a part of the electron diffraction profile, while the correspond- 
ing insets show a bigger range of momentum transfer with the prominent Bragg reflections as well 
as the superstructures of the basic crystallographic structure of 2i/-TaSe2 at low temperatures. 
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Single crystals of 2if-TaSe2 were grown from Ta metal and Se (99.95% and 99.999% 
purity, respectively) by iodine vapor transport in a gradient of 600 to 620 °C, the crystals 
growing in the cooler end of the sealed quartz tubes. A very slight excess of Se was included 
(typically 0.2 % of the charge) to ensure stoichiometry in the resulting crystals. Each exper- 
imental run lasted for 250 to 300 h. This procedure yielded single crystals with maximum 
size of 10 x 10 x 0.2 mm 3 . Appropriate thin films (thickness of about 100 nm) were prepared 
either by the use of an ultramicrotome or by cleaving with adhesive tape, whereby the crys- 
tal was repeatedly cleaved and the glue layer dissolved afterwards. The films were mounted 
onto standard electron microscopy grids and transferred to the spectrometer. As depicted 
in Fig.[TJ a very good sample quality is achieved, since discrete Bragg reflections as well as 
clear superstructures of the CDW arising at the predicted transition temperature (120 K) 
are visible. 

Potassium intercalation was achieved by evaporation from commercial SAES (SAES 
GETTERS S.P. A., Italy) getter sources under ultra high vacuum conditions (base pres- 
sure below 10~ 10 mbar). The intercalation levels were achieved by several discrete inter- 
calation steps, of which each took about 3 min. The relation between intercalation time 
and corresponding increase in potassium content could not be shown to be linear due to 
several reasons, such as varying getter quality and inhomogeneous saturation of the sample. 
Furthermore, the saturation potassium content was found to be stable over a time span of 
several months at room temperature and ultra-high vacuum conditions. All measurements 
were carried out with an 172 kV electron energy- loss spectrometer equipped with a He flow 
cryostat. The resolution of the spectrometer is 0.03 A -1 and 80meV for momentum and 
energy-loss, respectively. The spectrometer setup can be found in Ref. [26]. For the present 
studies, the access to elastic scattering processes by measuring at zero energy-loss is of great 
importance. 

To gauge and analyse the experimental results we performed density functional theory 
(DFT) calculations within the local-density-approximation (LDA). We used the scalar rela- 
tivistic mode of the full potential local orbital code (FPLO)[2Z] in version 9.09. In order to 
model the intercalation we used two approaches. In the virtual crystal approximation (VCA) 
K is+x-layers are inserted between adjacent TaSe2 layers, where 18 + x denotes the number 
of protons and electrons on the VCA-pottasium atom, which models the VCA equivalent of 
K^TaSe^ In this way the filling effect of the TaSe bands upon doping can be described. For 
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Figure 2. Loss- function at constant momentum transfer q showing a red shift of the 1 eV plasma 
resonance for increasing potassium content (left panel) . The right panel depicts the plasmon energy 
as well as the width of the plasmon peak, corresponding to the left panel, with respect to potassium 
content x per TaSe2 unit as estimated with the help of LDA calculations (cf. Fig.pp). 



judging the quality of VCA and for incorporating more realistic Fermi surface reconstruction 
effects we constructed several supercells (SC) with a corresponding K/Ta ratio. Guided by 
the experimentally observed superstructure patterns we focus on simple supercells with at 
least a threefold axis. For x = 1 /2 we also considered an orthorhombic cell. Especially for 
small concentrations it is not obvious that pottasium intercalates evenly into all interlayers. 
To cover such effects two series of supercells were created, one with potassium distributed 
as evenly as possible and one, where only every second interlayer contains the dopant (for 
x = ys, Y3, 1 /2). Details about the superstructures are found in the supplementary material. 
We did not try to determine the true structural ground state, since this could become rather 
involved and the purpose of the supercells is merely to supplement and check the VCA cal- 
culations. For all VCA and SC calculations we relaxed all lattice and Wyckoff parameters 
to accomodate breething of the lattice due to K-intercalation. The fc-integration for the self 
consistent calculations was performed via the linear tetrahedron method with a 12 x 12 x 4 
mesh for the simple VCA K2(i- x )Ta2Se4 unit cell. For the supercell calculations the num- 
ber of mesh points was adjusted accordingly to obtain an equivalent mesh density. Optical 
properties were calculated with a 32 x 32 x 24 mesh within the whole Brillouin zone, which 
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was checked for convergence. 

For the undoped compound we performed full-relativistic calculations to check the in- 
fluence of spin orbit coupling on the plasmon frequencies. It turns out that the spin-orbit- 
splitting of the bands at the Fermi level is quite sizable (400 meV) but the splittings (which 
reduce the z-dispersion) occur in such a way that their effects on the xy-plane plasmon ener- 
gies average out over the Brillouin zone. The resulting changes of in-plane u p are negligible, 
while the out-of-plane plasmon energy gets reduced. We only focus on the in-plane plasmon 
in this work and hence used scalar relativistic calculations throughout this work. 



RESULTS AND DISCUSSION 

Figure|2] (left panel) shows loss spectra of 2if-K x TaSe2 for increasing potassium content x. 
All spectra are normalized in the region around 2.5 eV and are furthermore corrected by fit- 
ting and subtracting the elastic line, contributing at around zero energy-loss [28]. The 
undoped sample shows a dominating peak at around leV that, with increasing potassium 
content, shifts to lower energies. Furthermore, the normalized peak maximum as well as 
the peak width are reduced. This spectral feature represents the charge carrier plasmon of 
the investigated samples, and therefore contains information on the charge carrier density. 
The reduction of the plasmon energy is a direct consequence of the charge carrier density 
variation induced by the potassium doping. Undoped 2if-TaSe2 is characterized by half 
filled conduction bands, and addition of either electrons or holes will lower the plasmon 
energy as has been discussed previously based on rigid band considerations [29]. In general, 
the plasmon energy therefore can be also used to analyze the doping level. Since the band 
structure of 2if-TaSe2 is rather complex and variations beyond a rigid band approach can- 
not be excluded, we compared calculated screened plasmon energies (Fig.|5j below) to the 
measured ones and deduced the x-axis in the right panel of Fig. [2} We used the thus defined 
concentrations in the rest of the paper. Note, however, that the calculated and measured 
screened plasmon energies are remarkably similar, which makes such a gauging procedure 
feasible (see also discussion below). 

It can be concluded that saturation of potassium doping is likely to be reached at about 
2i/-Ko.77TaSe2. Older reports on intercalation of transition-metal dichalcogenides were able 
to find also higher dopant rates up to x — 1. [131 EH] However, Rouxel [H] also emphasized 
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that a homogeneous structure in the selenides could only be obtained for rates of x = 0.6 
to 0.7. Furthermore, the spectral width of charge carrier plasmons in most cases is predom- 
inantly determined by damping via interband transitions [501 I3T] . As a consequence, we 
attribute the reduction of the plasmon width upon doping mainly to its shift to lower ener- 
gies, where less interband transitions are available for a damping of the plasmon excitation. 
In addition, mid-infrared transitions as observed for undoped 2iJ-TaSe2 are suppressed upon 
doping [32], which could also reduce the plasmon damping. 




Momentum Transfer (1 / A) 

Figure 3. Elastic scattering measurements of the afc-plane with increasing potassium content from 
top left to bottom right showing CDW superstructures for low temperatures and low K contents 
and arising dopant superstructures for increasing K contents. 

In Fig. [3] intensity maps, achieved by elastic scattering in EELS, are shown. They show 
congruent sections of the a6-plane with increasing potassium content from top (x = 0) to 
bottom (x = 0.77) at a temperature of 20 K. The bright spots in the highlighted TM and 
TK direction ([100] and [110], respectively) represent the basic structure of the compound. 
Additional reflections at one-third and two-third of the [100] vector are observable in the 
undoped phase (topmost left part of Fig.|3|, which can be assigned to the occurring CDW 
at temperatures below 120 K (see also Fig.[T]). [3H7j 

During the successive doping process several changes in the scattering image can be 
observed. The most prominent superstructure is a hexagon like ring at around 50 % doping 
(bottom left part of Fig.[3j). Upon increasing potassium addition, in this state the sample 
first shows a mixed superstructure of the CDW spots at one-third and two-third as well 
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Figure 4. Elastic scattering measurements of the fully doped sample (2f/~-Ko.77TaSe2) around the 
[110] zero reflection (left panel) showing the suppressed intensity for increasing angles of the scat- 
tering direction with respect to the a6-plane. The right panel shows the according [11/] reflections 
(see text), that emerge with increasing angles, allowing a calculation of the c-axis length, which is 
widened compared to the undoped crystal structure. 

as new superstructures occuring at one-half (for x = 0.15) and one and three quarter (for 
x = 0.64), respectively. However, the latter doping state is no longer showing the original 
triple CDW superstructure. At the highest possible doping level (for x = 0.77) again a ring 
like structure around the [100] reflection can be found. 

These findings lead to the conclusion that in 2if-TaSe2 the CDW is suppressed by K 
intercalation, as it was also found for the related compound 2if-TaS2 [33H35] . At the same 
time new features arise that must be assigned to an increasing potassium filling of the crystal 
structure. At this point a comparison between alkali metal intercalated graphite and TMDCs 
shows striking similarities. Both material classes show a van-der-Waals coupled hexagonal 
layer structure and metallic behavior in their room temperature phase |36| . Furthermore, a 
transition to a superconducting phase could be induced or supported, respectively, by means 
of alkali metal doping [35| 137] . Upon intercalation, most commonly Rb or K, graphite is 
able to form distinct phases, which can be distinguished as stage 1 and 2 (CsX and C14X, 
respectively, with X being the alkaline metal). With increasing content of the dopant it 
is showing a diversity of superstructures that can be assigned to the intercalant [38H4T] . 
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However, the formation of distinct phases present in intercalated graphite could not be 
observed in 2if-TaSe2. In particular the continuous variations in line shape and position of 
the charge carrier plasmon do not signal the simultaneous presence of two or more phases at 
any doping level. Nevertheless, the appearance of new superstructure peaks at one-half of 
the basic unit length of 2if-TaSe2 would be consistent with the potassium arrangement for 
a doping rate of x = Y 3 as used in our calculations (see supplementary informations Fig. 2 
left). Our theoretically assumed structures (see supplementary information) would also be 
in agreement with these diffraction peaks to persist up to relatively high doping levels. At 
very high doping the doping induced superstructure would then lead to diffraction peaks at 
one-third and two-third of the original positions in agreement with Fig.[3] (lower right panel). 
A quantitative analysis of our diffraction data is not possible because of multiple scattering 
effects. 

Since the momentum is transfered within the sample plane it is not possible in the 
present case to directly access the c-direction of the crystal. With the possibility to rotate the 
sample around an axis perpendicular to the beam direction, measurements are possible in a 
range of 55° to 90° with respect to this axis. Figure|4] shows elastic scattering measurements 
(fully doped sample) in a plane rotated relative to the sample plane by the denoted angle. 
While the left panel shows a clear decrease of the [110]-peak intensity (around 3.6 A -1 ) 
with increasing angles, the right panel shows rising features at higher momentum transfers, 
corresponding to reflections with a finite I component, such as 3.9 A -1 representing [114], 
4A _1 — [115] and 4.2A -1 — [116]. The corresponding I component of around 0.37 A -1 is in 
agreement with the value of an earlier X-ray investigation of doped 2if-TaSe2 [20], where 
this value was found for a maximum doping rate of x = 2 /z, which is in very good agreement 
to our stoichiometry analysis (see below). 

In order to quantitatively determine substantial structural changes, the behavior of the 
lattice parameters in the LDA calculation can be summarized by noting that the pla- 
nar lattice constant stays rather constant upon doping, while the c-axis expands with K- 
intercalation. Figure [5]b (black solid symbols) shows that in the supercell calculations there 
is a rather sudden widening of the lattice even on small amount of added potasium, due 
to the rigid ion cores, which need some space and due to the rigidity of the covalently 
bound TaSe2 prisms. For the cells, where only every second interlayer has K the c-axis is 
smaller (light (red) filled symbols), which follows from the same reasoning. With increasing 
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Figure 5. Theoretical results of the LDA calculations (virtual crystal approximation (VCA) as 
well as constructed supercells (SC)): a) £oo from the interband contributions; b) the c-axis lattice 
parameter; c) the plasmon peak position of the calculated loss function, the solid line shows a fit 
of the supercell frequencies, used for fixing the potassium content of the measured sample; d) the 
bare plasmon frequency, es denotes cells with alternatingly empty interlayers. 

K-content the c-axis shrinks because of the increasing electrostatic interaction provided by 
an increasing electron content in the interlayer. The behavior of VCA and SC coincide for 
larger concentrations. However, towards smaller x the c lattice parameter of VCA (open 
black and light(red) symbols) rises above the SC result, which can be understood from the 
fact that VCA contains the Kig-ion cores at a 100% filling rate for all concentrations, which 
keeps the TaSe2 layers appart, while for supercells all K-atoms vanish for x — > 0. This larger 
c-axis spacing however seems to mostly influence the plasma frequency along the z-direction, 
which is not considered in this investigation. However, the c-axis value of around 16 A for 
x ~ 0.8 and with it the corresponding I value (0.39 A -1 ) are slightly higher than the ones 
measured and from the literature (20], respectively. 

The planar bare plasmon frequencies of VCA and SC calculations Fig.[5]d form a consis- 
tent trend to lower energy values for increasing potassium content. We emphasize that for 
undoped 2if-TaSe2 the agreement between our measurements and the calculations is very 
good. Interestingly, the plasmon frequencies jump, when going from undoped to slightly 
doped samples, which we attribute to the immediate lattice widening upon small amounts 
of intercalated potassium. 
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In order to gain some insight into the screening due to interband transitions and with it 
the link to the measured plasma frequencies we also calculated the real part of the dielectric 
function in dependence on x. In general the real part tends to flatten out into a quasi 
constant (£«>) for w < leV. Unfortunately, there are some oscillations for small energies 
due to a low energy feature in Ime, whose strength is concentration dependent. This makes 
a clear extraction of Eoo rather approximate. Fig.[5^L shows the trend of Eoo as a function 
of x. The screening is of the same order as extracted from experiment and also exhibits 
the sudden jump, upon the doping-onset. Since £ ro is not a constant and since the real 
part has low energy features we used the theoretically calculated loss function to extract 
a screened plasmon peak position (Fig. [5^) and to estimate the K-content of the samples 
via the measured plasmons peak position. The values obtained in this manner are used 
throughout this contribution. 



SUMMARY 



To summarize, we performed Electron Energy-Loss Spectroscopy as well as LDA calcula- 
tions on the transition-metal dichalcogenide 2if-TaSe2 and thereby investigated the influence 
of in situ intercalation of potassium into the sample. We demonstrate that our intercalation 
procedure can be used to achieve rather high doping levels, which allows to study the physi- 
cal properties of these structures. Furthermore, we present a detailed theoretical analysis of 
the structure and electronic properties of intercalated compounds. Upon doping the charge 
carrier plasmon continuously shifts to lower energies, a behavior supported by the theoret- 
ical investigations and ruling out a possible formation of doped phases. The bands at the 
Fermi level, which get filled by potassium doping are of Ta5<i character. Their bandwidth 
W is rather small (below 2eV), which could hint towards the presence of sizable correlation 
effects. However, the Wannier functions for these band states as calculated in Ref.H2] are 
rather extended, which results in a comparatively small value for the Hubbard U. Hence, 
we expect a (possibly much) larger ratio w /u compared to 3d transistion metal compounds 
and no pronounced Mott-Hubbard physics. 

Structurally, our thin film samples show a transition to a charge-density wave state below 
120 K in the undoped case and a suppression of the according superstructure upon potas- 
sium intercalation. Moreover, new superstructures with a different periodicity arise that can 
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be assigned to the increasing potassium content in the samples. The changes in the basic 
structure with respect to the c-axis can be consistently found in the Electron Energy-Loss 
Spectroscopy measurements as well as in the calculations. The spectral width of the charge 
carrier plasmon significantly decreases upon potassium addition, i.e., the plasmon excita- 
tion is much less scattered at higher doping levels. Furthermore, as our data indicate the 
suppression of the charge density wave ground state upon electron doping, it is tempting to 
speculate about the possibility of an increasing superconducting transition temperature in 
potassium doped 2_£/-TaSe 2 (130 mK in the undoped compound [3]) following the observa- 
tion that often there is competition between a superconducting and a charge density wave 
state in transition metal dichalcogenides [3U EHl H3] . 
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SUPPLEMENTARY MATERIAL 

In order to check and complement the virtual crystal approximation (VCA) calculations 
we constructed a few supercells. No attempt was made to actually determine precise mini- 
mum energy structure for each concentration, which can get rather involved due to possible 
disorder effects or large wavelength structural modulations. We picked cells, which seemed 
reasonable and have a threefold axis (except for x = 1 /2, see below). For the smallest cell 
(50 % K) we tried three structures in order to get a feeling for the spread of the plasmon 
parameters and uJ Pbare due to structural and Fermi surface reconstruction effects. For 
the same reason we additionally considered cells with every second interlayer being empty 
for x = For x = l /& we only have an alternatingly empty supercell, due to cell size 
considerations. It turns out that the energies of the alternatingly empty cells are higher 
but only by a rather small amount of about 100 to 200 eV per TaSe2 formula unit. Such 
structures could be realized and stable in nature due to kinetic effects during intercalation 
and the large barrier for K to move through a TaSe2 layer. 

• • • c • c 
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Figure 1. Top and side view of the supercells for x = 0, i/s, Atoms at the corner of the prisms 
are Se, the atoms inside the prisms are Ta and the largest atoms on top are K. The difference 
between^ = i/s and x = 1 /a is the alternation of empty and filled K-layers for x = 1 /s. 
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Figure 2. Supercells for x = !/3. The cell with every second K-layer being empty is higher in enegry 
by 73mev per TaSe2 formula unit. 




Figure 3. Supercells for x = !/2. Note, that the second supercell has two different alternating 
K-layers. The first cell has orthorhombic symmetry. The first two cellsare energetically degenerate 
with a difference below 3 meV, while the energy of the alternatingly empty third cell os higher by 
193 meV per TaSe2 formula unit. 
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Figure 4. Supercells for x = 2 /3, 3 /4, 7 /8. Note, that for x = 7 /s there are two different alternating 
K-layers. For x = 1 (not shown) each TaSe2 prism is covered with one K. 
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